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Abstract 

The DENIS survey is currently imaging 21334 deg 2 of the mainly southern sky in the 
UK and the observations are expected to go on until mid 2000. The expectations for 
extragalactic and cosmological research are outlined, including a quantitative assessment 
of the effects of recent star formation on the measured fluxes of galaxies. The galaxy 
extraction is much improved with the modeling of the PSF across the 12' x 12' frames and 
the reliability of star /galaxy separation (currently based upon a combination of classical 
and neural-network based methods) is measured from visual inspection to be > 90% at 
I = 16. The I band counts follow the high bright-end normalization and the J differential 
counts follow N(J) ~ 11 x dex[0.6 (J — 14)] deg -2 mag -1 and are expected to be complete, 
reliable and photometrically accurate (Am < 0.1) for J < 14. 



1 Cosmology and Near-Infrared selection 

Workers in the Infrared bands know that the Near-Infrared (NIR) is 4 (J) to 10 (K) times 
less affected by extinction from dust than is the V band |J. This has two advantages: 1) A 
nearly-full view of the Universe, even behind the Zone of Avoidance, and 2) A view of external 
galaxies unaffected by their internal extinction by dust. 

The second advantage of the Near-Infrared is its smaller bias to recent star formation, in 
comparison with bluer bands and also to the mid and far IR, whose emission arises from thermal 
dust created and heated by young stars. 

The lack of sensitivity to recent star formation can be quantified as follows. Figure [I] shows 
the evolution of broad-band luminosities, normalized to that of a 13 Gyr old stellar population 
for different color bands. The evolution is strongest in the B band (and even more so in 
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Figure 1: Time evolution of broad band luminosities (Johnson BVJHK and Cousins RI) of 
stellar burst for stellar population with a Rana & Basu |f26| IMF, using the PEGASE [JL4J] 
spectral evolution model including nebular emission (which only affects the luminosities within 
10 Myr from the burst). 

U or far-UV, neither of which are plotted in Figure [I]) and weakest in the H and K bands. 
Nevertheless, there is a burst of emission in H and K roughly 10 Myr after the starburst, when 
the massive stars evolve off the main-sequence to the giant branch. Table [I] below shows that 
the enhancement of broad band flux in the B band due to star formation occurring 1 Myr ago, 
relative to the B flux of an old stellar population of the same mass and IMF, is 40 times greater 
than the analogous enhancement of the K band flux. 



Table 1: Broad-band luminosity enhancement relative to old (13 Gyr) population (normalized 
to the analogous K band relative enhancement) 



Age 


B 


V 


R c 


h 


J 


H 


K 


1 Myr 


40 


18 


14 


3.4 


1.8 


0.7 


1 


10 Myr 


5.6 


2.4 


1.7 


1.4 


1.06 


1.05 


1 


100 Myr 


7.8 


3.5 


2.3 


1.6 


1.15 


1.03 


1 


1 Gyr 


3.7 


2.7 


2.4 


1.9 


1.26 


1.04 


1 


2 Gyr 


2.1 


1.7 


1.5 


1.4 


1.11 


1.02 


1 



One interesting point is that the H band is to within a few percent as insensitive to recent 
star formation as the K band, and even less sensitive by up to 30% for very recent (< 5 Myr) 
star formation, with the inclusion of nebular emission in the adopted PEGASE |14| spectral 
evolution model, which is stronger in K because of free-free recombination enhancing the con- 
tinuum at longer wavelengths. Another interesting point is that the J band at 1.25 \i is almost 
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as insensitive to recent star formation as the K band, to within typically 20% or less, except for 
bursts younger than 8 Myr. The I c band is considerably more sensitive to recent star formation 
than J (in Figure p] and Table [I], Ic is typically three times more distant from K , logarithmically, 
than is J from K). 



2 The DENIS survey 

In December 1995, the DENIS (DEep Near Infrared Southern Sky Survey |]12|) began imaging 
the southern sky, delimited by —88° < 5 < +2°, in the I c (0.8 fi), J, and K s (iC-short at 
2.15 /i) wavebands. Its characteristics are given in Table ^| below. The data processing is done 
at the IAP (Paris) and Sterrewacht Leiden (Netherlands). The American 2MASS project [ffij 



will have a similar scope, with slightly more sensitive J and considerably more sensitive K 
channels, and an H channel instead of our / channel. 



Table 2: Characteristics of the DENIS survey (adapted from [|12 



Institutions 


Paris- Meudon, IAP, Leiden, Grenoble, Frascati, IAC, 
Innsbruck, Vienna, Lyon, Heidelberg, Budapest, 
Montpellier, Besangon, Sao Paolo 


Declinations 


-88° < 5 < +2° 


Telescope 


ESO lm (fully dedicated) 


Color bands 


I(0.8/im), J (1.25 H. ^(2.15/im) 

Tektronix CCD 1024 2 (I), Rockwell NICMOS-3 256 2 (JK) 


Detectors 


Pixel size 


I" (I), 3" (JK dithered at 1") 


Quantum efficiency 


0.65 (I), 0.8 (J), 0.61 (K) 


Exposure time 


9s (/), 8.8 x Is (JK) 


Read-out noise 


7e~ (I), 40 e~ (JK) 


Observing mode 


Stop & stare, concurrent 3-channel with dichroics 


Scan geometry 


12' (RA) x30° (Dec) 


Limiting magnitude (pt source 5a) 


18.0 (/), 16.0 (J), 13.5 (K) 


Saturation magnitude 


10.0 (I), 8.0 (J), 6.5 (K) 


Survey period 


December 1995 - mid 2000 


Data 


4000 GBytes (primary) 


Cost 


$3 million 



3 Extragalactic research and cosmology 

A variety of scientific questions concerning extragalactic science and cosmology will be answered 
with the catalogs coming out of the DENIS and 2MASS databases: 

Statistics of NIR properties of galaxies: DENIS and 2MASS will provide the first very 
large galaxy databases with NIR photometry, which should help for distance estimates of 
spiral galaxies with known velocity width (see ||3"2"|). Surface photometry of the brighter 



galaxies will be correlated with color as well as optical morphology, and with rotation 
velocity of known spirals |32| and used for distance estimates for ellipticals. 
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Cross-identification with other wavelengths: There will of course be plenty of cross-ident- 
ification of DENIS and 2MASS galaxies with samples at other wavelengths, such as optical 



galaxy samples (see |32|, [29]]), IRAS galaxies, quasars, radio-galaxies, galaxies found in 
blind HI surveys, etc. The NIR properties (mainly their location in color-color diagrams) 
of such objects will be targeted for discovering new large samples of such objects. There 
will naturally be followups at non-NIR wavelengths of DENIS and 2MASS galaxies (see 



Galaxy counts: From its J-band galaxy counts, DENIS should settle, once and for all, the 



debate on the bright-end of the galaxy counts, where first estimates |2T[] saw a depletion 
relative to the extrapolation of the faint-end counts, while later work || disputed this. 
At issue are the questions of galaxy evolution and whether the environment of the Local 
Group is underdense on very large scales [z ^ 0.1). 

Zone of avoidance: Studying galaxies behind the Galactic Plane has two main applications 
(see |^9|]): 1) Mapping the large-scale distribution of galaxies in this still poorly known 
region. Indeed, the Zone of Avoidance contains interesting structures such as the largest 
large-scale concentration of matter in the local Universe, the Great Attractor (at the 
intersection of the Supergalactic Plane and the Galactic Plane Jl6|) and within the Great 
Attractor, the Norma cluster, Abell 3627, richer and closer than the Coma cluster 117 



2) The fluxes and angular sizes of galaxies are affected by extinction from dust in the 
Galactic Plane, and one can measure this extinction from galaxy counts ||, colors [213 
and color- color diagrams |29[ 



Structures of galaxies: Only a few catalogs of clusters [T^, [13]] and compact groups [[25 



are based upon automatically selected galaxy samples, which happen to be optical and 
photographic (hence subject to photometric non-linearities). Because star formation is 
probably enhanced by galaxy interactions, one expects that the statistical properties of 
pairs, groups and clusters of galaxies built from NIR selected galaxy catalogs will be 
different from those built from optical catalogs. DENIS and 2MASS will thus have the 
double advantage of using an NIR galaxy sampled based upon linear (non-photographic) 
photometry. The applications of such NIR-based samples of structures of galaxies are 
and include understanding the dynamics of these structures, their bias to 



numerous 



projection effects, their constraints on Q and the primordial density fluctuation spectrum, 
their use as distance indicators, and the environmental influences on galaxies. 

The large-scale structure of the Universe: The NIR selection and the linear photometry 
will also benefit the measurement of statistics (two-point and higher-order angular corre- 
lation functions, counts in cells, topological genus, etc.) of the large-scale distribution of 
galaxies in the Universe. For example, the (3D) primordial density fluctuation spectrum 
of galaxy clustering can be obtained from the two-point angular correlation function || 
or from the 2D power spectrum ||. Moreover, by the end of DENIS and 2MASS, large- 
scale cosmological simulations with gas dynamics incorporated (thanks to which galaxies 
are properly identified) will provide adequate galaxy statistics in projection that will be 
compared with those obtained from the surveys, iterating over the cosmological input 
parameters of the simulations. 

DENIS has already made contributions in stellar astronomy, with a new short scale length 
(2.3 kpc) and cutoff (15 kpc) of the Galactic disk pTJ and the discovery of three candidate 



brown dwarfs [11], one of which (< 0.065 M in Corvus) has been confirmed spectroscopically 

m. 
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4 Galaxy extraction 

The advent of first DENIS results in stellar astronomy highlights the difficulty of extracting 
galaxies from the DENIS images. Indeed, at the star/galaxy separation limits of the / band 
(which is the most sensitive, see § §), there are 3 to 5 times more stars than galaxies at high 
galactic latitudes |T9[], and up to 10 4 times more in the Galactic Plane (with considerably 



less favorable ratios at brighter magnitudes, since galaxy counts rise faster than star counts). 
Whereas much of the science outlined in section § |^ requires galaxy samples that are complete 
and reliable (with few stars or optical defects misidentified as galaxies), such complete and 
reliable galaxy extraction is hindered by the current weaknesses in the DENIS camera: 

• Low sensitivity in K caused by the thermal emission of the instrument 

• Undersampled PSF in J and K (despite the dithering mentioned in Table 0) 

• Large PSF in / (typically 2"), caused in part by defocusing, coma and astigmatism of the 
camera 

• Important PSF variations across the frames (10% to 50%) with the same causes plus 
optical misalignment 

• Elongated PSFs near the frame edges 

• Artifacts caused by readout electronics 

In particular, the star/galaxy separation in the / band depends crucially on the accurate 
modeling of the PSF across the DENIS frames. An r.m.s. error of 5% on the FWHM of 
the PSF leads to a loss of 0.5 magnitude in the completeness limit of galaxy extraction at given 
reliability of star/galaxy separation |J. 

4.1 Pipeline 

We have worked on the images that were bias subtracted, flat-fielded, background subtracted, 
and bad pixel flagged with the software written by J. Borsenberger J7j, who has also performed 
absolute astrometry by tying to a copy at the CDS in Strasbourg of the USNO-A1.0 @ astro- 



metric catalog obtained from PMM scans. One of us (G.A.M.) has implemented a preliminary 
galaxy pipeline with the following steps: 

• Cosmic ray and bad pixel filtering 

• Reading photometric zero-points and airmasses from relevant files 

• Galaxy extraction using SExtractor version 1.2b6a (including a neural- network based 
star/galaxy separation ||, whose input parameters are 8 isophotal areas, the maximum 
intensity and as a control parameter, the FWHM of the PSF), with detection and Kron 
photometry [rj| parameters optimized from simulated images |I[ 



• Modeling of the variations of the object FWHMs across the frame (with a program written 
by M.T.) 

• Star/galaxy separation using the updated PSF at the object position obtained from the 
modeling across the frame, using E. Bertin's SExPlay2 

4.2 Photometric accuracy 

Photometry is an important aspect in understanding the diagnostics of galaxy extraction and of 
most algorithms of classical (not based on neural networks) star/galaxy separation. We report 
here on the first measures of photometric accuracy of DENIS candidate galaxies, using objects 
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lying in the 2' overlaps of adjacent images. Figure § shows the differences in photometry on 
overlap objects that are likely to be galaxies. The photometric accuracy on a given individual 




Figure 2: Photometric difference on overlap (concordant astrometry < 2") candidate galaxies 
on 50deg 2 , discounting flagged (blends, edge) objects and objects with centers within 20 pixels 
from edge of frames. The limiting neural network || stellarity index (0 for extended objects, 1 
for stars) for galaxies is given at the top (PSF modeling was used, see § fO| ). 

measurement is (mi — m2) rm s/2 1//2 . We thus infer from Figure ^] that 07 = 0.08 at / = 17, 
a j = 0.08 at J = 14.5, and a K = 0.42 (0.31) at K = 12 (12.5). These do not include systematic 
photometric errors caused by uncertain zero-points, extrapolations to total magnitude, etc. The 
small number of overlaps in K means that the photometric accuracy in K is uncertain. 



4.3 Star galaxy separation 

The optical deficiencies of the DENIS camera cause variations of the PSF across the DENIS 
frames, which results in an excess of objects classified as extended (galaxies) in regions where 
the PSF is larger, if no modeling of the PSF is performed before star/galaxy separation (23| . 



One of us (M.T.) has modeled the PSF FWHM as a two-dimensional second order polynomial 
across the frame using unflagged objects with / < 16.2 and iteratively rejecting 2.5 a outliers 
(see also M). 

Figure shows how the neural network stellarity seems to discriminate better between 
stars and galaxies once PSF modeling is incorporated, providing a gain of roughly one magni- 
tude in the reliability limit (the DENIS strip 5544 used here spans galactic latitudes from 15° 
to 34° towards the Galactic Center). 

Figure ^ shows the positions of candidate galaxies in the I band before and after modeling 
the PSF variations. Vignetting is apparent in both strips and corrected by the PSF modeling, 
perhaps too much so. Also, the upper left quadrant appears to produce an excess of galaxy 
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Figure 3: Neural network |5j stellarity parameter (1 = stellar, = extended) versus magnitude 
for 180 frames of fairly low latitude strip 5544, without (left) and with (right) PSF modeling. 



Figure 4: Positions of candidate galaxies (I < 16.0 and limiting neural network stellarity index 
for galaxies Sj < 0.1, only blending flags permitted) on DENIS frames, stacked over DENIS 
strips 5544 (top) and 5570 (bottom, 33° < b < 61°), each containing 180 images. The left and 
right plots correspond to neural network stellarities defined without and with PSF modeling, 
respectively. Circle diameters vary as the square root of the fluxes. 
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candidates before PSF modeling, which may also be present after PSF modeling. However, 
the positions of the extended objects appear to lack uniform coverage of the frame in the PSF 
modeled star/galaxy separation. Hence, the PSF modeling is still imperfect and can probably 
be improved. We are currently analyzing new ways of modeling the PSF for adequate uniformity 
of the extracted galaxies. 

One of us (W.B.) has obtained estimates of the reliability of star galaxy separation through 
systematic visual inspection of 56 objects in one 6deg 2 strip (5544), with / < 16.0 and who 
figure in the galaxy branch of the area-magnitude diagram shown in Figure |5|. Borderline cases 
were checked visually by a second person (G.A.M.). 91% of objects passing both classical area 
vs. magnitude and neural network star/galaxy separation were confirmed as galaxies, while 7% 
were double stars (half of these were in an unusually dense region) and 2% (1 object) was an 
optical default (the spike of a mildly saturated star). 



Figure 5: Area versus / magnitude for DENIS strip 5544 (objects at least 20 pixels from frame 
edges). Circles are candidate galaxies with stellarity Si < 0.10, the line is a fiducial classical 
star/galaxy separation, and the plus signs above the line are classified as stars by the neural 
network (5/ > 0.10). 

Similarly, we visually checked the 31 objects classified as stars by the neural network but 
situated above the classical dividing line in the area-magnitude diagram, and all were confirmed 
stars though 2 (6%) seemed to have faint galaxies superposed. For example, the isophotal area 
of the bright object at I = 12 classified as a star was seriously overestimated because it lied 
near a saturated star, but this did not fool the neural network. 
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5 Preliminary results on 50 deg z 



There are two reasons to rely on 7 band star/galaxy separation when building lists of galaxies 
in 7 or K: 1) 7 is considerably more sensitive than 7 and K (except in highly extinguished 
6 |p9|), and 2) 7 has better angular resolution. Since typical I — 7 colors 



regions, 



.4 



v 



> 



are almost always bluer than 1.6 (see below), 7 band star/galaxy separation will be roughly 
90% reliable (see § above) down to 7 = 16 — 1.6 = 14.4 and for K ^ 13.2 (for similar 
reasons), well beyond the completeness limits for galaxy detection (see Table below). 

The results below are based upon I band star/galaxy separation using stellarity (Sj < 0.10) 
and area vs. magnitude (dividing line in Figure |^). Figure ^| presents the 7,7 — 7 color- 
magnitude diagram for 8 high galactic latitude DENIS strips. The increased scatter at faint 7 
magnitude is caused by photometric errors. The very blue and very red objects are probably 
misclassified stars. In any event, only the points above the two lines in Figure |6| have reliable 
star/galaxy separation. The galaxy locus seems to occur in the region 0.6 < 7 — 7 < 1.6. At 
7 = 14.0 (14.4), 96% (90%) of the galaxies have 7 < 16.0, so that this 7 list is over 87% (82%) 
reliable. 
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Figure 6: Color magnitude diagram for galaxies from 8 high latitude DENIS strips (Si < 0.1 and 
isophotoal area in / above critical line of Figure ||, with no additional star/galaxy separation in 
J; objects closer than 50 pixels to the frame edges are excluded). The horizontal line represents 
7 = 16 and the tilted line represents J = 14.5. 



Note that the magnitudes in these plots are based upon Kron [18] photometry, which at- 
tempts to approximate the total magnitude of galaxies. Since, 7 is more sensitive than 7, the 
Kron aperture radius is larger in 7. We thus expect fixed aperture colors to appear slightly 
redder. 

In figure |7|, we show the UK color-color diagram for our 8 high latitude strips. The 
galaxies with reliable photometry [large filled circles in Figure [?D cluster around 7 — 7 ~ 1.2 
and 7 — K ~ 1.1, and an important part of the scatter is caused by photometric errors. Again, 
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Figure 7: Color-color diagram for galaxies (using the same star/galaxy separation as in Figure 
^ and no additional star/galaxy separation in K) in 8 high latitude DENIS strips. The large 
circles correspond to J < 16.0 and K < 11.5, the small circles to I < 16.0 and K > 11.5, the 
dots to I > 16. The solid and open circles correspond to J < 14.4 and J > 14.4, respectively. 




Figure 8: DENIS differential galaxy counts (using the same star/galaxy separation as in Figure 
|, but with an additional color criterion deduced from Figure U 0.6 < I — J < 1.6). Objects 
closer than 50 pixels to the frame edges are excluded, for a total equivalent area of 50 deg 2 . The 
dotted (left), dashed (middle), and solid (right) histograms correspond to the DENIS counts in 
K, J, and I c , respectively. Data points are counts from Gardner et al. [15j in K (squares) and 
I c (filled circles) and from Lidman & Peterson [|19[] in I c (open circles). 
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the colors within fixed apertures will be redder (see above). We are in the process of visualizing 
the galaxies that appear unusually red in J — K, but we note that most of these are too faint 
for very reliable star/galaxy separation in / and/or accurate photometry in J or K. 

Figure || shows the galaxy counts in UK bands for the 8 high galactic latitude DENIS 
strips. The DENIS / counts follow Gardner et a/.'s [15] out to a 75% completeness limit of 
/ ~ 17.25. The disagreement at the bright end with Lidman & Peterson [19|] argues for a high 
bright-end normalization of galaxy counts (see ||), although not as high as Gardner et a/.'s. 

Note that DENIS, Gardner et al. and Lidman & Peterson all work with the Cousins / 
band, so no conversion was made from another I filter. Also, our survey has smaller error 
bars at the bright end as it covers 4 to 5 times the solid angle of the two cited surveys. In 
comparison with Gardner et a/.'s [fO^ K counts, the DENIS K counts become incomplete at 
K ~ 11.25. Figure [8] shows what we believe are the first reliable bright galaxy counts in J, and 
are roughly 80% (50%) complete out to J ~ 15.25 (15.75). The J counts can be parameterized 
by N(J) ~ 11 x dex[0.6 (J - 14)] deg^mag^ 1 . 



6 Discussion 

With the plots presented above, we are now in a position to provide our expected limits for the 
DENIS survey. From Table ^| below, we infer that star/galaxy separation is the limiting factor 
in / and J, while completeness and accurate (0.20 magnitude) photometry are the limiting 
factors in K. And from the counts of Figure |], we estimate that complete, reliable, and 



Table 3: Estimated DENIS limits from 50 deg 2 of reduced data 





h 


J 


K 


completeness (~ 80%) 


17.25 


15.25 


11.25 


star/galaxy separation (90% reliability, from J) 


16.0 


14.0 


12.5 


photometry (0.10 magnitude accuracy) 


17.4 


15.0 


<11? 



photometrically accurate DENIS catalogs will include 600 000, 200 000, and < 8000 galaxies 
at / < 16.0, J < 14.0, and K < 11.2, respectively. A cooling system has been very recently 
installed on the DENIS camera, and we hope that this will improve the i^-band sensitivity. 
We also expect that better star/galaxy separation, in particular through better modeling of the 
PSF variations across the frames, will increase the reliability limit up to / ~ 16.5, which will 
roughly double the number of galaxies in the complete and reliable / and J band catalogs. 
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